Matching rules from Al-Co potentials in an almost realistic model 
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We consider a model decagonal quasicrystal of composition Also.iCoig.g - closely related to actual struc- 
tures, and using realistic pair potentials - on a quasilattice of candidate sites. Its ground state, according to 
simulations, is a Hexagon-Boat-Star tiling satisfying Penrose's matching rules. In this note, we rationalize these 
results in terms of the potentials; the Al-Co second-neighbor potential well is crucial. 



From the discovery of quasicrystals, it has been difficult 
to decide whether they are stabilized by entropy or by en- 
ergy. (Here "energy" means an ideally quasiperiodic ground 
state.) In the decagonal case, the fastest route to a resolu- 
tion may be ab-initio-based modeling. Recently, we observed 
a ground state [JJ that perfectly implements Penrose's match- 
ing rules [1^ in a toy model of a binary Al-Co quasicrystal, 
closely related to realistic d(AlNiCo) models ||l,|j,|5t]. In this 
contribution, we sketch the matching rules' origin; in particu- 
lar, how fine-tuned must the potentials be to obtain this result? 

Our simulations use a recipe introduced to study real AI- 
Co-Ni phases iH 0, HI The only input data are the number 
density, the composition ratio, and the (quasi)lattice constants 
(tile edge is = 2.455 A and layer spacing c/2 = 2.04A).In 
an initial "unconstrained" Monte Carlo simulation the atoms 
hop as a lattice gas on discrete, properly placed candidate 
sites (see Pii4]), which decorate random rhombus tilings (that 
are rearranged as another kind of Monte Carlo move). The 
lowest energy configuration visited during the run is saved 
(which is far better than a typical state, due to energy fluc- 
tuations in these relatively small systems). The second stage 
is a "constrained" simulation flipping HBS tiles with a fixed 
atomic decoration (which allows larger systems) inferred from 
the unconstrained results. 

The result (Fig. [1) is a Hexagon-Boat-Star (HBS) tiling of 
edge 2.45 A with Al atoms at exterior vertices and Co atoms 
on the interior rhombus vertex, as well as (respectively) one, 
two, or two internal Al atoms. Importantly, an Al vertex site 
goes vacant if it would have no Co nearest neighbors. For the 
present work, we introduced a "half-constrained" simulation 
in which atoms hop as a lattice gas on HBS tiles, limited to a 
site list of the ideal vertex sites, plus sites out along the 
midline of each Fat rhombus. This revealed that rearrange- 
ments are induced next to the empty vertices, such that the 
ground state is not literally the Penrose tiling (e.g. the de- 
fect marked by a pair of V's in Fig. [T]) However, we believe 
quasiperiodicity is maintained (the shape of hyperatoms in 5D 
space will change slightly). 
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FIG. 1: Best unconstrained configuration from original discovery 
runs, placing atom content AI169C042 in a smallish cell 31.9 x 
23.3 A, with periodic boundary conditions. We cooled gradually 
to T « HOOK with 10^ trial atom swaps per swappable pair, and 
2000 tile flips per flippable pair during the run. Co and Al atoms 
are shown by black and gray circles; large/small circles indicate the 
top/bottom layers. The seven V-rule violations (i.e. fat rhombus/thin 
rhombus) are marked "V", the three fat/fat violations are marked with 
circles, and the three vacated Al sites are marked by asterisks. Four 
violations are necessitated by periodic boundary condition. The pair 
of V's in the lower right corner mark an energetically favored defect 
(see text). 



I. ORDEIUNG AND DIAGNOSTICS 

The interactions (we used "GPT" potentials IH |^) have 
the following properties: (i) Vaiai(^) has a strong hard- 
core repulsion at close neighbor distances (R < 2.8A). (ii) 
Vaico(^) is extremely attractive at i? « 2.5A; (iii) VcoCo(^) 
has a rather strong attraction at the second-neighbor distance 
R « 4. 5 A, due to the prominent Friedel oscillations; (iv) sim- 
ilarly Vaico(^) is attractive, too, around R w 4.5A; this well 
is mainly responsible for the matching rules 

How does this order follow from these potentials? The 
question falls into two parts, (i) Why the HBS tiling? (ii) 
Granted the atoms do form an HBS tiling, why do the interac- 
tions force matching rules? This is the focus of Sec. [Ill which 
uses inflated (edge 4 A) HBS tiles to examine the energies 
with more detail than Ref . 1 1 ] . 

In Ref. yj we identified four aspects of ordering, the first 
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TABLE I: Important bond distances and potential values, bonds of 
each kind N{R) per cell, and change AE^j, in the contribution to the 
total energy change E{Ti,) — E{Ta), per cell. (Interlayer distances 
are written R2.) The same cell and atom composition were used as 
in Fig. [T] temperatures were « 5000,K Ti fa 2000K (using a 
randomly chosen typical configuration after equilibration, from the 
"half-constrained" simlation); To refers to the near ground state, the 
best found during a run at low temperature. 



R 


pair 


Vab{R) 


N{R) 


-ASo.i 




(A) 


{AB) 


(eV) 




(eV/cell) 


(eV/cell) 


2.455 


Al-Al 


0.4124 


119 


-1-0.8247 


-1-5.7729 




Al-Co 


-0.2919 


152 


-0.5839 


-1-1.4597 


2.542 


Al-Al 


0.2758 


248 


-1.1033 


-3.8615 




Al-Co 


-0.2583 


240 


-F2.0661 


-1-1.0331 




Co-Co 


0.1001 





0.0000 


-1-0.2002 


2.886 


Al-Al 


0.0791 


170 


-0.7119 


-0.7119 




Al-Co 


-0.0427 


52 


-^0.0427 


-1-0.3416 




Co-Co 


0.0950 





-^0.0000 


-1-0.3801 


3.016 


Al-Al 


0.0669 





-fO.0669 


-1-0.0000 


3.192 


Al-Al 


0.0637 


274 


-fO.6373 


-1-1.5296 




Al-Co 


0.0913 





-1-0.9127 


-1-0.0000 


3.789 


Al-Al 


0.0264 


126 


-0.0528 


-0.5281 




Al-Co 


0.0554 


2 


-fO.5538 


-1-1.4953 


3.846 


Al-Al 


0.0212 


126 


-fO.6375 


-1-0.3825 




Al-Co 


0.0424 





-fO.0848 


-1-1.2722 


4.303 


Al-Al 


-0.0055 


144 


-0.0439 


-0.0439 


4.465 


Al-Al 


-0.0071 


456 


-^0.3811 


-1-0.4940 




Al-Co 


-0.0349 


528 


-Fl.3261 


-1-1.3959 




Co-Co 


-0.0909 


116 


-0.5453 


-1-1.2724 


4.669 


Al-Al 


-0.0049 


283 


-fO.0147 


-0.0978 




Al-Co 


-0.0260 


149 


-^0.0521 


-1-0.8337 




Co-Co 


-0.0796 


45 


-^0.3 183 


-0.1591 


4.762 


Al-Al 


-0.0028 


247 


-0.1039 


-0.2254 




Al-Co 


-0.0189 


304 


-0.0756 


-1-0.1702 


4.997 


Al-Al 


0.0015 


340 


-0.0278 


-0.0278 




Al-Co 


0.0012 


104 


-0.0023 


-0.0187 




Co-Co 


-0.0193 





0.0000 


-0.1543 



two of which develop at relatively high temperatures (T > 
Ti = 1000 K). Aspect (1) is formation of the 2.45-A HBS 
tiling, with Co in the center and Al decorating the corners of 
every tile, as is quite universal in Al-transition metal decago- 
nals. HBS tiles trivially satisfy the Penrose double arrows. 
The HBS formation is attributed to the nearest-neighbor Al- 
Al and Al-Co interactions |1]. Medium-temperature config- 
urations from an "unconstrained" simulation may have many 
defects (of matching rules, or placement of Al atoms internal 
to the HBS tiles), but the HBS framework itself is inviolate. 

Aspect (2) is the Co-Co network due to second-neighbor 
Co-Co attraction: a supertiling of inflated HBS tiles (edge 
raij « 4 A) and small generalizations (e.g. pillow tile, be- 
low). But at medium T, the Co-Co tiling is not yet a Penrose 
inflation of the 2.45A-level HBS tiling. 

The other two aspects of order appear at T < Ti and im- 
plement matching rules. Aspect (3) is the "V-rule": a convex 
(27r/5) corner of the Thin rhombus in the (2. 45 A) H or B tile 
can only adjoin a concave (Svr/S) corner of the B or S tile, and 
vice versa. This ties together the HBS tiling and the Co super- 
tiling, requiring that every Co-Co edge have a Fat Hexagon 
around it. The V-rule is simply the Penrose (single-arrow) 



matching rule between Fat and Thin rhombi, as it is mani- 
fested on our (small) HBS tiles. Aspect (4) of the ordering - 
the final stage - is to satisfy Penrose single arrows between 
two Fat rhombi, on the edges not involved in the V-rule (two 
on every Hexagon and Boat); its explanation is deferred to 

Sec. In] 

Of course, the V-rule demands an equal number of the re- 
spective kinds of corner, thus ny — 2n/f + ub — 2nB + 
5ns- Combined with the formula for the numbers of rhombi, 
uh + Sue + 5ns = ntat and 2nH + ub = ntmn, this 
completely specifies the number of H, B, and S tiles. In the 
limit of zero phason strain, we get densities {nn ,nB,ns) = 
(r~^, r~^, r~^/\/5) (per rhombus) exactly as in the Pen- 
rose tiUng. (Also, if matching rules are obeyed, the vacated 
Al sites are centers of inflated Star tiles and their density is 
nvac = T~'^ / \/5, which gives the Also.iCoig.g stoichiometry, 
and the point density of 0.0697/A^.) 

There were two key diagnostics for our initial identification 
of the key interactions for matching rules. First, we compared 
the contributions of each pair distance R to the total energy 
as a function of temperature (Table IJl. As evidence for the 
story just given, big changes are seen upon cooling Too Ti, 
but none from Ti downwards, in (i) nearest neighbor [2.45- 
2.86A] Al-Al and Al-Co, forming HBS tiles, and (ii) in Co- 
Co 4.47^.67 A, forming the Co-Co network. The big changes 
at lower T were in Al-Co at 3. 79-4. 67 A, in accord with the 
analysis of Sec. HIl (below). 

The second diagnostic was to vary the potential cutoff ra- 
dius Tcut (which was usually 7 A). At j'cut — 5.lA the Pen- 
rose rules are still satisfied whereas at Tcut ~ 3.5 A one 
gets an HBS tiling but no matching rules, confirming indepen- 
dently that HBS formation is due to nearest-neighbor bonds 
but matching rules come from second-neighbor bonds. Fur- 
thermore^jf Tcut = 5.lA and internal Al's are removed from 
all tiles, pTOf] the tiling implements the V-rule but not the Fat- 
Fat rule, confirming the role assigned in Sec. to the "6" 
atoms in Fig. |2] and similar Al. 



II. MATCfflNG RULES VIA 4.0 A TILING 

Aspect (4) of the ordering - the final stage - is to satisfy 
Penrose arrows on the remaining Fat rhombus edges (two on 
every Hexagon and Boat). Let's focus now on the Fat rhom- 
bus/Fat rhombus matching rule, the crucial aspect (4) of the 
ordering process. To check the energy cost of a Fat/Fat viola- 
tion, which depends on the internal Al in those Fat rhombi. 
We reformulate the matching-rule problem in terms of the 
4.0A supertiles, in which arrows are defined on every edge 
by the Al atom position. A given supertiling can be broken 
into HBS tiles in several ways, corresponding to different ar- 
rowings of its edges. Note that usually (as in our 2.45A-HBS 
tiling), the degrees of freedom are in the tiles, and possible 
mismatches occur along tile-tile edges; whereas in the 4.0A- 
HBS model, each edge is an independent degree of freedom, 
and the possible mismatches occur between two edges in the 
same tile. (The same was true for the near-matching-rule of 
d(AlCuCo) 101.) 
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(a). (b). (c). 




FIG. 2: Inflated edges (broad shaded lines), as forced by the "V- 
rule". Sohd hnes show (a) H, (b) B, and (c) S tile; dashed lines 
are fragments forced by the V-rule. Insets below show the arrow 
decorations at the inflated scale (d). The arrowings of edges marked 
"?" are undetermined by the V-rule. Interactions between the labeled 
atom sites determine the Penrose arrowing as optimal, (e). Besides 
H, B, and S, only this pillow tile is consistent with the "V-rule", but 
it forces the unfavorable Al-Al distance of 3.79A. 



The Co-Co network has 4.46A bonds, along supertile 
edges, and 4.67A bonds. All of the latter relate endpoints 
of a 27r/5 supertile corner (e.g. in Fig.|2c)). Assuming the 
V-rule, [see Fig.|2a,b,c)], that makes five bonds per (small) B 
and five per small S tile. Hence, if the H-B-S content is fixed 
(as implied by the V-rule), the number of 4.67A Co-Co bonds 
is fixed and they do not contribute to the matching rule. (That 
same number, in the super tiling, is 2 per H, 3 per B, and 5 per 
S tile, so the super HBS content is also constrained.) 

What supertiles are possible? A corollary of the V-rule is 
that, at any (27r/5) [resp. 3(27r/5)] corner, both aiTows point 
out from [resp. into] the corner. The other constraint is that no 
interior space be left after a supertile's border is decorated by 
HBS tiles (filling an interior always violates the V-rule). The 
conclusion is we can have the super-H, B, and S tiles, or the 
"pillow" tile shown in Fig.|2|e) (or extensions of it by adding 
more segments of alternating 2(27r/5) and 3(27r/5) corners.) 

So let's consider the super-Boat tile (Fig. |2]i (d)). The V- 
rule leaves undetermined whether the two edges marked "?" 
both point to the same corner (as shown) or are flipped with 
one in and one out (this version violates the Penrose rule). De- 
pending on the choice for the first arrow, we either have Al at 
ai and 6i, or else at a2 and 62; similarly the primed sites de- 
pend on the second arrow. The energy difference between the 
options for first aiTow depends on 16 possible interactions be- 
tween sites fli, 61, 02, ^2 and the fixed atoms marked 1,2,3,4. 
(Several other fixed atoms interact, but they are symmetrically 
placed relative to the flip, so they can't affect the energy differ- 



ence.) In addition, there are four interactions between {ai, hi) 
and {a'i, 6^, however we must exclude those that connect two 
sites on the same tile. (That term was counted aheady in the 
energy per tile, and as noted earlier the total counts of H, B, 
and S tiles were fixed by the V-rule.) 

Here are the biggest contributions (energies are in Table 
either option makes one Al-Co(4.67) bond; the "right" 
way also gets two Al-Co(4.47) - quite favorable - and one 
A1-A1(2.89) - very unfavorable; the "wrong" way gets four 
A1-A1(3.85) and two A1-A1(3.79) which are unfavorable. The 
"right" way also makes eight A1-A1(4.47) - favorable - ver- 
sus just two for the "wrong" way. The net difference is 0. 1 825 
eV in favor of the "right" way. Preliminary numerical tests, 
correlating total energy with the number of matching rule vio- 
lations, indicate a typical matching-rule defect (of either kind) 
actually costs the order of 0.1 eV. 

Next, consider the super Hexagon. Its side edge arrows 
are free to point either way, but a calculation like that for the 
Boat says these arrows are favored to be parallel, which is the 
proper Penrose arrowing of a Boat. Finally, the V-rule permits 
just one possible arrowing on the super-Pillow tile (Fig. |2l), 
which forces a Fat-Fat rhombus violation in the middle asso- 
ciated with unfavorable Al-Al distances. Hence we expect this 
tile loses out to the super-Star (which carries the same count 
of2.45A HBS tiles). 



III. DISCUSSION 

The matching rules are mostly implemented by Al-Co inter- 
action at ^ 3.8 A and 4.462 A; however, they are a resultant 
of many terms, and the cost of a rule violation seems to be 
context-dependent (especially with our standard potential cut- 
off at 7A, which includes the third well.) Fortunately, as seen 
in the example of Sec. [Ill this isn't a "frustrated" problem: 
most contributions have the same signs. We checked the cost 
£^match of a rule violation: a fit of total energy versus number 
of violations (counted by hand) in 10 low-energy configura- 
tions (from T w lO'^K) gave i^match ~ O.leV. If we took a 
(near) ground state and moved the violation along a "worm" 
by local HBS flips, however, the cost was ^ 0.5eV, suggest- 
ing the minimum four violations forced by periodic boundary 
conditions had found sites where their cost was anomalously 
low. 

The matching rule depends on having the exact ratio of 
H:B:S tiles. That is a worry: by local flips, one can trade tiles 
HS BB or vice versa. Let's define energies EH^ Eb, Es 
as the sum of interactions within each tile (including a chem- 
ical potential /iai, since a flip changes the number of Al 
atoms). Then one naively expects to maximize or minimize 
ub, depending whether Eh + Es — 2Eb is positive or neg- 
ative, which would prevent reaching the right composition. 
But the tile-tile (matching-rule) interactions create a "gap" in 
Al site energies so that running the reaction in either direction 
costs a positive energy. 

The rule also depends on having the right internal Al dec- 
oration in every HBS tile. We studied the same model with 
Al-Ni: those potentials are quite similar to Al-Co (except the 
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Al-Ni attraction is not quite as strong as Al-Co). The Penrose 
tiling may well be the ground state for Al-Ni — indeed, the 
constrained simulation behaves the same for Al-Co and Al- 
Ni — but it was less robust in the unconstrained simulation: 
e.g., B tiles with just one internal Al on the mirror axis are 
metastable. 

As a reality check, we note that according to pair poten- 
tials, the (relaxed) matching-rule structure is unstable by 69.9 
meV/atom compared to the tie-line between AlgCo2 (19% Co) 
and other competing phases. [For this calculation we used 
the optimal "half-constrained" simulation result for the same 
size as Fig.[T]] The (relaxed) ab-initio energy, using the VASP 
package [8], came out unstable by 98.5 meV/atom. 1 1 1] 

For comparison, the best Al-Co structure models are 
« 10 meV/atom above the tie-line. In typical decagonal ap- 
proximants, the energy is reduced ^ 50meV/atom after Al 
atoms are allowed to "pucker" out of the layers (doubling the 
c axis to 8A), and to organize the proper correlations between 
the puckerings of nearby atoms. We have not yet investi- 
gated puckering in the present system, which would require 
molecular-dynamics simulation followed by relaxation ||j,|5t]. 



We also attempted to make our structure more realistic by 
going to a ternary, having recognized that the internal sites in 
HBS tiles (besides the interior vertex) are always "problem 
sites" in Al-transition metal decagonals. (E.g., those are the 
atoms that "pucker" |i5J.) The root problem is that the two Al 
in one Boat are a bit overpacked, in view of the Al-Al hardcore 
distance. Could we fill these sites with (somewhat smaller) 
Cu atoms? No: when we tried an Al-Cu-Co ternary (with pair 
potentials), Cu atoms entered Hexagons forming CuCo pairs 
(as seen earlier in Ref. 
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